DESY 09-025 
Edinburgh 2009/02 
LTH 822 



Electroproduction of the A^*(1535) resonance at large momentum transfer 
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We report on the first lattice calculation of light-cone distribution amplitudes of the Af*(1535) resonance, 
which are used to calculate the transition form factors at large momentum transfers using light-cone sum rules. 
In the region > 2 GeV^, where the light-cone expansion is expected to converge, the results appear to be in 
good agreement with the experimental data. 



Introduction. — Electroproduction of nucleon resonances 
has long been recognized as an important tool in the explo- 
ration of the nucleon structure at different scales. Transition 
form factors to nucleon excited states allow one to study the 
relevant degrees of freedom, the wave function and the inter- 
action between the constituents. Quantum chromodynamics 
(QCD) predicts [1, 2, 3, 4, 5, 6] that at large momentum trans- 
fer the form factors become increasingly dominated by the 
contribution of the valence Fock state with small transverse 
separation between the partons. There is a growing consensus 
that perturbative QCD (pQCD) factorization based on hard 
gluon exchange is not reached at present energies; however, 
the emergence of quarks and gluons as the adequate degrees 
of freedom is expected to happen earlier, at ~ a few GeV^. 
In this transition region the form factors can be measured to 
high accuracy, see e.g. [ /], and such measurements are in fact 
part of the experimental proposal for the 12 GeV upgrade at 
Jefferson Lab [8]. 

Theoretical progress in the QCD description of transition 
form factors has been slow. The major problem is that any at- 
tempt at a quantitative description of form factors in the tran- 
sition region must include soft nonperturbative contributions 
which correspond to the overlap integrals of the soft wave 
functions, see, e.g., [9, 10]. In particular, models of general- 
ized parton distributions (GPDs) usually are chosen such that 
the experimental data on form factors are described by the soft 
contributions alone, cf. [11, 12, 13]. A subtle point for these 
semi-phenomenological approaches is to avoid double count- 
ing of hard rescattering contributions "hidden" in the model- 
dependent hadron wave functions or GPD parametrizations. 
An approach that is more directly connected to QCD is based 
on the light-cone sum rules (LCSRs) [14, 15]. This technique 
is attractive because in LCSRs "soft" contributions to the form 
factors are calculated as an expansion in terms of the momen- 
tum fraction distributions of partons at small transverse sepa- 
rations, dubbed distribution amplitudes (DAs), which are the 
same quantities that enter the calculation in pQCD, and there 



is no double counting. Thus the LCSRs provide one with the 
most direct relation of the hadron form factors and DAs that is 
available at present, with no other nonperturbative parameters. 
Unfortunately, with the exception of the A(1232) resonance, 
up to now there exists almost no information on the DAs of 
nucleon resonances. Thus pQCD predictions [16, 17] cannot 
be quantified and the LCSRs cannot be used as well. 

Moments of the DAs can, however, be calculated on the lat- 
tice. In this work we suggest a synthetic approach combining 
the constraints on DAs from a lattice calculation with LCSRs 
to calculate the form factors. As the first demonstration of 
this strategy we consider the electroproduction of iV*(1535), 
the parity partner of the nucleon. This is a special case be- 
cause lattice calculations of baryon correlation functions al- 
ways yield results for baryons of both parities, = 1/2+ 
and = 1/2^ (see, e.g., [18, 19]), so in fact the results 
for iV*(1535) appear to be a byproduct of our calculation of 
the nucleon DAs [20, 21], to which we refer for further tech- 
nical details. We find that the shapes of the nucleon and N* 
DAs are rather different. A preliminary account of this study 
was presented in | |. In this paper we further use our results 
on the DAs to calculate the helicity amplitudes Ai/2(Q^) and 
S'i/2((9^) for the electroproduction of iV*(1535) in the LCSR 
approach. In the region > 2 GeV^, where the light-cone 
expansion may be expected to converge, the results appear to 
be in good agreement with the experimental data. 

Distribution Amplitudes. — The leading-twist(=3) nu- 
cleon (proton) DA can be defined from a matrix element of a 
nonlocal light-ray operator that involves quark fields of given 
helicity gTli) = (1/2)(1 ±75)9 [23]: 

= -\fNP-n^ul{P)j [dx] e-'^-" ^ Vn{x,).{\) 

Here P^, P^ = mf^, is the proton momentum, un{P) 
the usual Dirac spinor in relativistic normalization, an 
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arbitrary light-like vector = and G the charge-con- 
jugation matrix. The variables 2:1,2:2, 2:3 have the mean- 
ing of the momentum fractions carried by the three valence 
quarks and the integration measure is defined as j[dx\ = 
dxidx2dx36{^ Xi — 1). The Wilson lines that ensure 
gauge invariance are inserted between the quarks; they are not 
shown for brevity. 

The nonlocal operator on the l.h.s. of (1) does not have 
a definite parity. Thus the same operator couples also to 
A^* (1535) and one can define the corresponding leading-twist 
DA as 

{0\e''''- (uj iain)CMj (02/1)) fidl (agn) |iV* (P)) 
= IfN'P-riflu^j,, (P) j [dx] e-'^-" ^ (x,) , 



where, of course, 

and /at. are defined as 



The normahzation constants 



ijk 



(u,C^Mj)(0)T5#dfe(0)|iV(P))=/jvP -n^UNiP) 
{0\e'^''{u,C^u,){Qyf5^dk{Q)\N*{P)) 

= fN'P-n-i^iiUN>{P). (2) 

On the lattice one can calculate moments of the DA 



dx] x\x^2^3 fi^i) : 



which are related to matrix elements of local three-quark op- 
erators with covariant derivatives, see [2 1 ] for details. The 
normalization is such that ip'^'^^ = 1. 

There exist three independent subleading twist-4 distribu- 
tion amplitudes $^ , ^'l^ , E,^ (as for the nucleon). They 
can be defined as (cf. [23, 24]) 



\e'^'[ul {ain)Ci^u'^ia2n)jfdlia3n)\N*iP)) 
= ^P-nf (P) j [dx] e"'^" ^ 



uliain)C^-f^Mjia2n)j-f^i^dlia3n)\N*{P)) 
= ^^P -nfimN'U^N.iP)! [dx] e-»-P-" 2: ^.a. 

X 'fN^^f'^^'^ix.O-Xl^fix,)] , 
(0 1 e'^'' (uj {ain)Cfiiu] (aan)) ^4 (agn) | N* (P)) 

= ^P-nfi mi,, (P) j[dx] e-*^" ^ Sf * (x.) , 



where ^f''^^{x^) and * 



N',WW I 



4 yxi) aiiu >i-4 {xi) are the so-called 
Wandzura-Wilczek contributions, which can be expressed in 
terms of the leading-twist DA [ ]. The two new normaliza- 
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3.234(63)(86) 
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35.57(65)(136) 
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0.1311(113)(382) 


0.1034(160)(584) 
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~ 0.095 


0.0613(89)(319) 


0.0398(132)(497) 
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21 


~ 0.095 


0.1091(41)(152) 


0.1281(56)(131) 




2 

21 


~ 0.095 


0.1092(67)(219) 


0.1210(101)(304) 



TABLE I: The normalization constants (in units of 10 ^ GeV^) and 
moments of the DAs obtained from QCDSF/DIK configurations at 
P = 5.40 for the nucleon (iV) and TV* (1535) at = 1 GeV^ 
The first error is statistical, the second error represents the uncer- 
tainty due to the chiral extrapolation and renormalization. The sys- 
tematic error should be considered with caution. 



tion constants are given by the local matrix elements 

{0[e'^''{u,Cj^u,){Oyf,rdkmN*{P)) 

= A*TOAr.75UAr.(P), 

(0|e*^'*'-("^Ca^.%-)(0)75^^"4(0)|A^*(P)) 

= XlniN'JbUN'iP)- 

The asymptotic distribution amplitudes (at very large 
scales) for the nucleon and N* are the same: 

0^"(xO = 120x1x2x3 , $4'(xO = 24xiX2 , 

$f^--(x,) = 24x1X2(1 + ^(1-5x3)), 

M/f '^■-(x,) = 24x1X3(1 + ^(1 - 5x2)) , 
S4(xj) = 24x2X3 , *|'(xj) = 24xiX3 . 

Baryon states of different parity can be identified in a lat- 
tice calculation as those propagating forward or backward 
in (imaginary) time as long as their momentum vanishes 
[18, 19]. While vanishing momentum is sufficient for the eval- 
uation of the normalization constants, the higher moments of 
the DAs require nonzero momentum. In this case the signal 
in the negative parity channel is contaminated by a contri- 
bution of the ~ 1/2+ (nucleon) ground state of the or- 
der / {mNmN-) enhanced by the factor e'^-^™*~^™)'^^~*\ 
where T is the time extent of our lattice. However, this ef- 
fect seems to be quite small in our results: Replacing the par- 
ity projector (1/2) (1 + 74) by (1/2) (1 + irriN/ENh^) [18] 
changes the first (second) moments of the N* DAs by 1% 
(5%), which is well below the statistical error. In principle, 
there is still a contamination by the = 1/2+ A^*(1440) 
(Roper) resonance, but for small momenta this effect is ex- 
pected to be negligible [18]. Another issue is that in the 
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physical spectrum there are two = l/2~ resonances, 
A^*(1535) and iV*(1650), which are close to each other, so 
that they cannot be distinguished by means of their mass dif- 
ference in our calculation. Because of the peculiar decay pat- 
tern of A^*(1650) we expect, however, that this state has a 
much smaller coupling to the usual interpolating operators 
[19]. So our results can be identified with the contribution 
of TV* (1535) alone. All these questions certainly deserve a 
further study. 

The results of our calculation of the normalization constants 
fw. A*, A2 and of the first few moments of the leading-twist 
DA of the A^*(1535) resonance are compared to the similar 
calculation for the nucleon ' ] in Table I. It attracts attention 
that /a?* is about 50% larger than /^r. This means that the 
wave function of the three quarks at the origin is larger in the 
= 1~ state than in the = 1+ state, which may be coun- 
terintuitive. The momentum fraction carried by the u-quark 
with the same helicity as the baryon itself, (p^^", appears to 
be considerably larger for N*, indicating that its DA is more 
asymmetric. These features suggest that the large asymme- 
try of the nucleon DA observed in QCD sum rule calculations 
[25, 26, 27] may be due to a contamination of the sum rules 
by states of opposite parity, which are difficult to separate in 
this approach. 

The calculated moments can be used to model the A^* 
leading-twist DA as an expansion in orthogonal polynomials 
corresponding to the contributions of multiplicatively renor- 
malizable operators (in leading order), see [ ]. The compari- 
son of such models for N and N* , obtained using the polyno- 
mial expansion to second order and the central values of the 
lattice parameters, is shown in Fig. 1. 

Helicity Amplitudes from LCSRs. — The matrix element 
of the electromagnetic current j"" between spin- 1/2 states of 
opposite parity can be parametrized in terms of two indepen- 
dent form factors, which can be chosen as 

{N*{P')\jT'\N{P)) = um.{P')i^T,un{P), 

r,. = — n — {(iq-y - q : — o-z^pr , (3) 



rriN 



where q = P' — P is the momentum transfer. The helicity 
amplitudes Ai/2((3^) and S'i/2(Q^) for the electroproduction 
of A^*(1535) can be expressed in terms of these form factors 
[28]: 



A, 



/2 



eB 

e 



^1/2 = -j=BC {mN~mN^)Gi{Q') + mNG2{Q' 
V 2 L 

Here e is the elementary charge and B, C are kinematic fac- 
tors defined as 
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FIG. 1: Barycentric plot of the distribution amplitudes for nucleon 
(up) and A^* (1535) (down) at i-ijfg = IGeV using the central values 
of the lattice results. The lines of constant xi, X2 and xs are parallel 
to the sides of the triangle labelled by X2, x^, and x\, respectively. 



The LCSRs are derived from the correlation function 



dxe'^''^N*{P)\T{lm3T{x)m, 



where 77 is a suitable operator with nucleon quantum numbers, 
e.g. the loffe current [ U]. Making use of the duality of QCD 
quark-gluon and hadronic degrees of freedom through disper- 
sion relations one can write a representation for the form fac- 
tors appearing in (3) in terms of the DAs of N*. In leading or- 
der, the sum rules for Q^Gi{Q^)/{mNmN') and — 2G'2(Q^) 
have the same functional form as the similar sum rules [14, 15] 
for the Dirac and Pauli electromagnetic form factors of the 
proton, with the replacement rriN rriM' in the light cone 
expansion part, and different DAs. 

In the present calculation we used a model for the leading- 
twist DA including first order corrections in the polynomial 
expansion, asymptotic expressions for the "genuine" twist- 
4 DAs and the corresponding Wandzura-Wilczek corrections 
up to twist-6 as given in Ref. [23]. The results are shown in 
Fig. 2. The shaded areas correspond to the uncertainty in the 
lattice values as given in Table I. In the region > 2 GeV^, 
where the light-cone expansion may be expected to converge, 
the results appear to be in good agreement with the experi- 
mental data. 
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FIG. 2; The LCSR calculation for the helicity amplitudes Ai/2{Q^) 
and 5'i/2(Q^) for the electroproduction of the N*{1535) resonance 
using the lattice results from Table I for the lowest moments of the 
A^*(1535) DAs. The curves are obtained using the central values of 
the lattice parameters, and the shaded areas show the corresponding 
uncertainty. 



Discussion and Conclusions. — In this work we suggest 
to calculate transition form factors for nucleon resonances at 
intermediate momentum tranfer, combining the constraints on 
DAs from a lattice calculation with LCSRs This approach 
seems to be especially promising for N*{1535), the parity 
partner of the nucleon, because of the relative ease to sepa- 
rate the states of opposite parity on the lattice. The accuracy 
is expected to increase significantly when calculations with 
smaller pion masses and on larger lattices become available. 
This would remove a major source of uncertainties which is 
due to the chiral extrapolation. 

In order to match the accuracy of the lattice results, the 
LCSR calculations of baryon form factors will have to be ad- 
vanced to include NLO radiative corrections, as it has become 
standard for meson decays. For the first effort in this direction 
see [35]. In addition, one needs a technique for the resumma- 
tion of "kinematic" corrections to the sum rules that are due 
to the large masses of the resonances. 
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